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I — I Abstract 

We study the 'tilting' of Mott insulators of bosons into metastable states. These are described 
by Hamiltonians acting on resonant subspaces, and have rich possibilities for correlated phases 
with non-trivial entanglement of pseudospin degrees of freedom encoded in the boson density. 
^ , We extend a previous study (Phys. Rev. B 66, 075128 (2002)) of cubic lattices to a variety of 

lattices and tilt directions in 2 dimensions: square, decorated square, triangular, and kagome. 
For certain configurations three-body interactions are necessary to ensure that the energy of the 
^ effective resonant subspace is bounded from below. We find quantum phases with Ising density 

0. wave order, with superfluidity transverse to the tilt direction, and a quantum liquid state with no 

o 

C_) broken symmetry. The existence of the quantum liquids state is shown by an exact solution for a 

particular correlated boson model. We also find cases for which the resonant subspace is described 
j> by effective quantum dimer models. 
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I. INTRODUCTION 



The observation of the superfluid-insulator quantum phase transition of ultracold bosonic 
atoms in an optical lattic^ has launched a wide effort to engineering other correlated quan- 
tum states of trapped atoms. Much of the effort has focused on entangling the spin quantum 
number of the atoms between different lattice sites: superexchange interactions between 
neighboring lattice sites have been observed^, but longer range spin correlations have not 
yet been achieved. 

Here we focus on a proposaP^lto obtain quantum correlated phases in a pseudospin degree 
of freedom which measures changes in boson number across the links of a tilted lattic^^. 
The interactions which entangle the pseudospins are not suppressed by factors oit/U (where 
t is the tunneling between lattice sites, and U is the local boson repulsion), and so non-trivial 
entanglement is likely to be readily accessibldP. This proposal takes advantage of the unique 
feature of ultracold atom systems, namely, the ability to study many-body non-equilibrium 
physic^^'^^'^, and suggests that various intriguing phases can appear as metastable stated. 
Moreover, the recent advent of site-resolved imaging of ultracold atoms in optical lattice^^^^ 
opens the way to the direct detection of the pseudospin and the investigation of its quantum 
correlations. 

The previous theoretical work on tilted lattice^^ focused on simple cubic lattices in 
spatial dimensions d = 1, 2, 3, with the tilt along one of the principal cubic axes. In d = 1, 
this theory predicted Ising density wave order for strong tilt, and this has been observed 
in recent experiment^. Here, we will study a variety of other lattices and tilt directions. 
We will find analogs of the density wave and anisotropic superfiuid phases found earlier. 
However we will also find new quantum liquid states, and some cases which map onto the 
quantum dimer model. These states entangle the displacements of all the atoms, and so 
should be attractive targets of future experiments. We will find an exact liquid ground state 
for a particular correlated boson model. 

We begin by recalling the Hamiltonian of a tilted Mott insulator. It is described by the 
generalized bosonic Hubbard modeP with an additional potential gradient along a certain 
direction: H = i^kin + -f^u + -f^tiit: 

H^n = -tYl ij>lh + %) (1-la) 

<ij> 

Hv = |5^ni(n,-l) + ^5^n,(r2,-l)(ni-2) + ... (1.1b) 

i i 

Htilt = -E'^e-riUi. (1.1c) 

i 

Here bi are canonical boson operators on lattice sites i at spatial co-ordinate r^, and hi = h\hi. 
The first term in H\] describes two body interactions. The second term is an effective three 



2 



o 



o 



o 



o 



(a) 



o 



o 



(b) 



(c) 



FIG. 1: (a) Parent Mott insulator, which is the ground state for A 
Ising density wave order for A — )• — oo. 



oo. (b,c) Ground states with 



body interaction, generated by virtual processes involving higher band^^^*^. Such a term is 
present in ultracold atomic systems, as has recently been measurecP^**^. As we will show, 
the presence of this term dramatically changes the physics of tilted two dimensional lattices, 
therefore we include it in our model. The potential gradient is E, and the fixed vector e is 
normalized so that the smallest change in potential energy between neighboring lattice sites 
has magnitude E. We will assume a filling of one atom per site in the parent Mott insulator 
in most of the following discussion. 

We are interested in the dynamics of the resonant subspace that appears when the lattice 
is tilted so that E becomes of order U. Specifically, we will work in the limit 



|A| ,t < E,U, 



[1.2) 



where we define 



We also define our tuning parameter 



A = U-E. 
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and allow A to range over all real values. Thus we include processes which carry energy 
denominators of order U — E to all orders, but neglect processes which have energy denom- 
inators of order U or E. 

Let us now review the properties^ of a tilted chain of sites in ci = 1; see Fig. [T| As 



explained in Ref. |3], the limit Eq. (1.2) defines a strongly interacting many-body problem 



in the resonant subspace, which can be described by an effective Hamiltonian. Unlike the 



underlying Hubbard model in Eq. (1.1c), the energy of this effective Hamiltonian is bounded 
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from below, and it makes sense to study its stable ground state and low-lying excitations 
by the traditional methods of equilibrium many body theory. This will help us describe the 
dynamics of experimental studies in the metastable resonant subspace of the full Hubbard 
model in Eq. (1.1c) defined by Eq. (1.2). It was found that the resonant subspace was 
described by a quantum Ising model for an Ising pseudospin residing on the links of the 
chain: an up spin represented the transfer of an atom across the link (or the creation 
of a dipole particle-hole pair), while the down spin had no change from the parent Mott 
insulator. In other words, in a parent Mott insulator with occupancy Uq per site, a down 
spin was a pair of sites with occupancy (no, no), while an up spin was (no — 1, no + 1). The 
Ising model had both a longitudinal and transverse external field, representing the energy 
of the dipole and the tunneling amplitude of the atom respectively. However, two nearest 
neighbor dipoles cannot be created simultaneously, and this implied the presence of infinitely 
strong nearest-neighbor antiferromagnetic exchange interactions. It is this infinitely strong 
exchange which allows easy access to strong-correlation physics in such tilted latticed. The 
phase diagram of the Ising antiferromagnet with both longitudinal and transverse fields has 
also been studied by others,'^^'^ with different physical motivations. For A — oo, the Ising 
model had a paramagnetic ground state, corresponding to configurations close to the parent 
Mott insulator. In contrast for A — )• — oo, the Ising model had antiferromagnetic long- 
range order, corresponding to a maximum density of particle-hole pair excitations about 
the Mott insulator. There are two possible ways of arranging these excitations, depending 
upon whether the particles reside on the odd or even sites, and this ordering is linked to the 
antiferromagnetic order of the Ising model. In between these limit quantum phase 

transition was founcP at an intermediate value of A: this transition was in the universality 
class of the quantum Ising ferromagnet in a transverse field and zero longitudinal fielcP'2IE21_ 
Next, let us revieW^ cubic lattices in dimensions d > 1, with a tilt along a principal lattice 
axis, with \Us\ ^ t. Then we have to distinguish the physics parallel and transverse to the 
tilt direction. Parallel to the tilt direction, the physics is similar to d = 1: there is an Ising 
density wave order associated with dipolar particle-hole pairs, which turns on as A decreases. 
In contrast, transverse to the tilt direction, the quantum motion of individual particles and 
holes is allowed, not just of dipole bound states. Such motion raises the possibility of 
Bose condensation of particles/holes and of the appearance of superfluidity. However, the 
single particle or hole motion is constrained to be purely transverse. Consequently there 
are an infinite number of conservation laws, one for each lattice layer orthogonal to the tilt 
direction: the total number of particles in any layer is constrained to equal the total number 
of holes in the preceding layer. Because of these conservation laws, the superfluidity is not 
global, but restricted to each transverse layer separately: there is no Josephson tunneling 
term which links the superfiuid order parameter of adjacent layers. Of course, because of 
the translational symmetry in the effective Hamiltonian parallel to the tilt, the superfluidity 
appears simultaneously in all the layers. In this superfiuid state, atomic motion is insulating 
parallel to the tilt, and superfiuid transverse to the tilt. Such transverse-superfluid states also 
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appear as A decreases, and can co-exist with the Ising order parallel to the tilt. The previous 
workP did not note the U3 term, even though its presence is required for the applicability of 
the d > 1 dimensional results therein. 

In the present paper we will examine several tilted two dimensional lattices. We find 
that the physics not only depends on the lattice geometry and tilt direction, but also on the 
magnitude of U3 compared to other energy scales of the systems. This U3 term is needed 
to guarantee the stability of some lattices, as for example the cubic lattices reviewed above. 
The reason is the following. Let us assume a parent Mott insulator with filling uq = 1. 
Processes of the kind (1, 1) — )■ (0, 2) cost repulsion energy U, and gain tilt-energy E (if 
the two participating sites differ in potential energy E). These processes are tuned to be 
resonant. Depending on the lattice structure, it may then happen that two sites differing 
in potential E are each occupied by two bosons. The process (2, 2) — t- (1, 3) costs repulsion 
energy U + f/3, and gains potential energy E. In the case {Usl ^ t this process is off- 
resonant. We can then identify a resonant subspace which focuses on the (1,1) — (0,2) 
processes only. This was implicitly assumed in Ref. [3] without discussion of the crucial 
role played by U3, and will be assumed in our Sections [ITj|IV[ In such cases, our effective 
HamiltoniarP^ method applies, the energy of the resonant subspace is bounded from below, 
and we can use equilibrium methods to study the phase diagrams. 

On the other hand, if lU^l < t, we cannot neglect processes of the type (2,2) — )■ (1,3). 
In the presence of such processes, the resonant Hilbert space is increased. This can have 
dramatic consequences on the physics; in some cases, we find that the energy in the resonant 
subspace becomes unbounded from below, and the system becomes unstable. 

The present paper will examine several lattice configurations of Eq. (1.1c) in c? = 2. We 
will assume jf/al ^ t in all cases, except in Section |Vj We will find analogs of the Ising and 
transverse-superfiuid phases described above in these lattices. However, we will also find 
a qualitatively new phenomenon: the appearance of quantum liquid ground states with no 
broken symmetries. These phases appear in these lattices for two main reasons. First, even 
for a large tilt, not all the sites can participate in the creation of particles or holes due to the 
geometrical constraints present in these lattice configurations. Such geometrical constraints 
generate exponentially large degenerate states in the limit of infinite tilt. Secondly, the 
lattice structure ensures that there is no free motion of unbound particles or holes in the 
direction transverse to the tilt in the resonant subspace. The absence of such unbound 
motion prevents transverse-superfluidity, allowing for the appearance of the quantum liquid 
states. 

The outline of the remainder of the paper is as follows. In Section |Tl] we study a square 
lattice with a tilt in a diagonal direction, while in Section III we study a triangular lattice 
with different tilt directions. In these cases, we will find analogs of phases studied previously, 
involving Ising density wave and transverse-superfiuid order. In Section IV we will consider 
the kagome and certain decorated square lattices. Here, the transverse-superfiuid order is 
suppressed for certain tilt directions, and novel quantum liquid ground states appear. Also, 



5 



FIG. 2: (Color online) Full blue circles denote quasiparticles, and black circles with light yellow 
filling denote quasiholes. The tilt vector is e = (1, 1) and the resonant particle-hole excitations are 
as shown in the figure. 

for another tilt direction on the kagome, the Ising density wave order is suppressed, and 
we obtain effective decoupled one dimensional systems. In Section |V| we briefly discuss the 
1/3 = case. 



II. SQUARE LATTICE 

As noted in Section [l| this section always assumes \Us\ S> t. We have already reviewed 
the results of Ref . E] for a square lattice tilted along a principle lattice direction. Ref . [3] also 
briefly considered a tilt along a diagonal lattice direction, and we will discuss this case more 
completely in the present section. 



For the diagonal tilt, we choose the vector e = (1, 1) in the Hubbard model in Eq. (1.1c). 
In this situation each site has two neighbors to which resonant tunneling is possible; see 
Fig. [2} Here and in the following, we label the lattice sites as (/,m) where / (m) represents 
the X (y) coordinate, respectively. Using the methods of Ref. [3l the Hamiltonian of the 
resonant subspace can be written as 



Hph = -t^yno{no + l)J2 [{p\m+i) + Pli+i)m) + Hx. +^Y1 {pLpim + h\Ji 



X / I I L I / 

''Im 

Im Lm 



(2.1) 

Here we have introduced bosonic quasiparticles p]^ and quasiholes h]^, and we identify the 
parent Mott insulator with flUing no, \Mno), as quasiparticle and quasihole vacuum, |0), 



and so 



hl^ |0) := -^km \Mno) . (2.3) 
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FIG. 3: (Color online) The matrix elements for the above processes are +teff for process (a), and 
— teff for (b) and (c). 



These operators are hard-core bosons and so we have the constraints 



VlmVl^yn < 1, 



hLhir. 



vlmVlrnni^riim 



0. 



(2.4a) 
(2.4b) 
(2.4c) 



We are now interested in describing the global ground state of the effective Hamiltonian in 



Eq. (2.1) while it is subject to the constraints in Eqs. (2.4). It is useful to first consider the 



limits of A — 7- oo and A — — oo, followed by a general discussion. 



A. Limit A — )• cxd 



To zeroth order in 1/A, the unique ground state is the vacuum, the parent Mott insulator. 
All particle-hole excitations are gapped with the same energy (see Fig. [2]), Ei = A. At second 
order in 1/A, we obtain an effective Hamiltonian for particles and holes, whose structure is 



strongly constrained by Eq. (2.4). A number of distinct processes appear in the perturbation 
theory: 

• Overall energy shift due to virtual processes. The vacuum couples to all states with 
one neighboring particle-hole pair, and so the vacuum energy is shifted down. Each 
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particle-hole state couples to zero and two particle-hole states. To second order the 
energy Eq for the vacuum and Ei for states with one particle hole pair are: 



Eo = -2iV-noK + l), (2.5) 



E, = =A-i2N- 8)^no(no + 1). (2.6) 

where N is the number of lattice sites. 

• Hopping of particles and holes along direction transverse to the tilt: Particles and 
holes can hop individually via second order processes along lines which are orthogonal 
to the tilt direction, see Fig. |3)3,c. The magnitude of the effective hopping is 

t^nojno + 1) 

|A| ■ ^ 

It is this process which leads to transverse superfluidity. Of course, because the su- 
perfluidity is now one- dimensional, it is only quasi-long range, and characteristic of a 
Luttinger liquid. 

• When particles and holes are proximate to each other, we should consider their hopping 
as contributing to the hopping of a dipole particle-hole pair. One such process is shown 
in Fig. [3^, and it leads to the motion of a particle (or rotation of a dipole) with matrix 
element opposite to that without an adjacent hole. 

Collecting these processes, we can write down an effective Hamiltonian in the manifold of 
excited states with one particle hole pair (with energy ^ A). We label {Ip, nip) the position 
of the quasiparticle, and (Z/^, nih) the position of the hole. Clearly the vacuum state only 
couples to states where particle and hole are on neighboring transverse diagonals, so that 
'^p + lp = + Ih + ^ = d. We only need three integers to describe the positions of particle 
and hole: d,mp,mh. We label the states by 



\d, TUp, rrih) = pLm„.m„^Lm,-l,m, |0) (2-^ 



-m,p,m,p d'-mii — lj'mf, 

and the effective Hamiltonian in this manifold is 



-f^ph, cs = El + tefT 1^ {\d, m, m) + \d, m,m — 1)) ((d + 1, m, m| + (ci + 1, m + 1, m\) 

d,m 

+ \d, m, m) {d, m + 1, m| + \d, m, m) {d,m,m — 1\ + H.c. 
-tes ^ ^\d,mp,mh) {d,mp + l,mh\ + \d,mp,mh) {d,mp,mh - l\ +B..C. . (2.9) 
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FIG. 4: (Color online) A — )• — oo: Three possible states with maximum dipole density. At zeroth 
order in 1/|A|, these states are degenerate, along with an infinite number of other state with 
the same dipole density. At second order in 1/|A|, two states are selected as the global ground 
states: the checkerboard state on the top-left, and its symmetry-related partner. The state in the 
top-right has a domain wall between the two checkerboard states: the domain wall costs energy 
t^/{2 |C/ - S|) per unit length. 

Here the first sum is for hopping of particle-hole pairs, both in transverse and longitudinal 
direction of the tilt. The second sum is for particles and holes hopping individually in 
transverse direction. Diagonalization of this Hamiltonian, as in Ref. 131 yields a continuum 
of separated particle-hole excitations, along with a dipolar particle-hole bound state. 

B. Limit A — > — oo 

At zeroth order in 1/|A| there are an infinite number of degenerate ground states, all 
maximizing the number of particle-hole pairs. Three such states are shown in Fig. |4j In 
Ref. 131 it was stated that these degenerate ground states can be labelled by the set of dimer 
coverings of the square lattice: this is incorrect, because given a set of particle and hole 
configurations, there is, in general, no unique assignment of them into nearest-neighbor 
dipoles. 

We can examine the lifting of the ground state degeneracy in a perturbation theory in 
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1/1 A|. At leading order, the matrix elements in the ground state subspace are all diagonal. 
Each particle that has one hole to the left or below, can undergo a virtual annihilation, 
reducing the energy by (t^/|A|)?T,o(?T-o + 1) = tcs- If it has two neighbors that are holes (left 
and below) it can do two virtual hopping processes, reducing its energy by 2teff. Thus we 
should optimize the latter configurations: this leads to two degenerate ground states, which 
look like checkerboards. These states break a lattice translation symmetry, and there is an 
Ising order parameter. 

C. Discussion 

The basic phenomenology that has emerged from our discussion so far of the diagonal 
tilted square lattice is quite similar to that for the principal axis tilt considered in Ref. |3l 
For A — )■ oo, we have continua of particle-hole excitations, while for A — t- — oo we have a state 
with Ising density wave order, and domain wall excitations. While for the principal axis tilt, 
these results appeared already at first order in 1/|A|, here for the diagonal tilt case we had 
to go to order 1/ A^ to obtain the free hopping of particles and holes, and the appearance of 
the Ising order. We don't expect these distinctions to be important at smaller values of |A|, 
and so the phase diagram of the diagonal tilt case should be similar to that of the principal 
axis tilt: Ising order appearing as A decreases, along with transverse quasi-superfluidity at 
intermediate values of A. 

III. TRIANGULAR LATTICE 

As in Section |ll| this section also assumes lU^] S> t. The triangular lattice case is 
very similar to the square lattice case. The phase transitions are to Ising density-wave 
ordered states and possibly transverse quasi-superfiuid phases. The tilt breaks the rotation 
symmetry, so the dipole states on the tilted triangular lattice are not frustrated. 

A. Tilt along a principal lattice direction 

For the tilt along a lattice direction, e.g. ai, we choose the tilt magnitude so that the 
creations of dipoles along the two other lattice directions, 02 and 03, are resonant (see Fig.|5^) 
but the dipole creation along the direction Si is not. Therefore, the effective Hamiltionian 
in the resonant subspace for this case is the same as for the square lattice with a diagonal 
tilt, and all conclusions from Section [TT] apply here. 

One could also choose the tilt magnitude so that processes along di are resonant. Then 
all processes along 02 and 03 are off-resonant, and the resonant subspace separates into 
decoupled one- dimensional systems. 
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(b) 

FIG. 5: (Color online) (a) Triangular lattice tilted along principal lattice direction. The effective 
Hamiltonian of the resonant subspace for this system is the same as the one for a square lattice with 
diagonal tilt. Full blue circles denote quasiparticles, and black circles with yellow filling denote 
quasiholes. (b) Triangular lattice tilted perpendicular to a principal lattice direction. For large 
negative A, there are two- fold degenerate ground states with broken lattice symmetry. 

B. Tilt perpendicular to a principal lattice direction 

Here we briefly consider a triangular lattice with tilt along a principal direction. In the 
limit A — J- oo, the ground state is the parent Mott insulator and all excitations are gapped. 
Particles and holes can separate and hop along the axis transverse to the tilt direction and 
reduce their kinetic energy. 

On the other hand, for large and negative A, there are precisely two degenerate ground 
states, as illustrated in Fig. |5}3. There is an Ising density wave order parameter associated 
with these states. These states and their symmetries are similar to those found for the tilted 
square lattice in Ref. |3l and a similar phase diagram is expected. 



IV. DECORATED SQUARE AND KAGOME LATTICES 

this section also assumes \Us\ ^ t. We will now consider 
frustrated models in the sense that in the strongly tilted limit, A — ?■ — oo, not all sites can 
participate in dipole creation. This requires lattices with a larger unit cell, as will become 
clear from our discussion. Furthermore, the tilt direction of those lattices can be chosen such 
that independent motion of particles and holes is not possible in the resonant subspace. We 
can then label the resonant subspace by a set of dipole coverings. In the limit A — > — oo, 
we will obtain a large degeneracy in the dense dipole states. This degeneracy is lifted by 
corrections in inverse powers of |A|. The leading corrections are not diagonal in the basis of 
dipole coverings: this sets up the possibility for novel quantum liquid phases. 



As in Sections O and III 



11 















> i 
1 — • — ( 


> 

» — 


c 
c 

> ( 
( — • — ( 


» 




J 


— ( 




) 

> 1 


> < 


» ( 















e 



(a) 

— ►e 

(b) 

FIG. 6: (Color online) (a) Decorated square lattice in diagonal tilt, and (b) kagome lattice with tilt 
along lattice direction, where the tilt strength is chosen so that only processes along the diagonal 
lines are resonant. These two models lead to the same effective Hamiltonian in the resonant 
subspace. Excitations of the parent Mott insulator are dipoles (shown as think lines connecting 
quasiparticle and quasihole). The particles and holes cannot separate. Each dipole link blocks four 
other links from forming dipoles (shown as crosses). Those are also the four links onto which each 
dipole is allowed to hop through virtual processes. 

A. Decorated square lattice with a diagonal tilt 

Here we consider a decorated square lattice tilted along a diagonal direction, as illustrated 
in Fig. [6^. We define the distance between nearest-neighbor sites to be a so that the unit 
cell of the square lattice in Fig. [6^ has size 2a. 

Interestingly, the effective Hamiltonian for the resonant subspace resulting from the tilted 
kagome lattice shown in Fig. [6]d is the same as the one obtained in the decorated square 
lattice in Fig. [6^. While the geometry of the kagome lattice is different from the decorated 
square lattice, for the specific tilt of Fig. |6]d, dipoles cannot be resonantly created along the 
links parallel to the tilt direction. This reduces the geometry of the tilted kagome to the 
tilted decorated square lattice. Using this identity, we will present most of the results using 



12 



the decorated square lattice illustrated in Fig. [6^. 

The decorated square lattice has 3 sites per unit cell. Let be the number of unit cells, 
and we assume periodic boundary conditions. There are two kinds of sites on this lattice: 

sites that have four neighbors, and 2N sites that only have two neighbors. We define the 
unit cell so that a site with four neighbors is in its center. We will refer to these sites as 
central sites. 

The Hamiltonian in the resonant subspace is 

H = ^Y^ dida - t^noK + 1) XI {'^'^ + ' 

a a 

where a labels the lattice links, and dl creates a dipole on this link. The resonant subspace 
has the constraints 

dida < 1, (4.1) 

and 

didadiJa' = (4.2) 

if a and a' share a lattice site, see Fig. |6] Note that particles and holes cannot separate 
on this lattice through hopping in the resonant subspace. It is notable that this effective 
model restores the rotational symmetry of the lattice and is symmetric under rotation of 
7r/2 around the central sites. 



1. Limit A — >■ oo 



Here the ground state is the parent Mott insulating state. The excitations are gapped 
and the lowest excitations are given by the creation of a single dipole. To zeroth order in 
1/A, all states with one dipole are degenerate. We now describe the effective Hamiltonian 
for the manifold of excited states with energy A to second order in 1/A. First of all, there 
are energy shifts for the vacuum energy Eq and energy of one dipole Ei, given by 

Eo = -4iVteff (4.3) 

El = A-{m- 6%s (4.4) 

with 

t^nojno + 1) 

= ^] • f^^-^) 

Additionally, second order processes lead to hopping of the dipoles. Note that a dipole 
always involves one central site and one neighboring site located to the right, left, up or 
down relative to the central site. Thus, we label the dipole states by \l,m,a), where (/,m) 
denotes the position of the central site and a = 1,2,3,4 indicates the direction of the dipole, 
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FIG. 7: (Color online) Labeling of sites and links on the decorated square lattice. 



see Fig. [7} With this notation, the effective Hamiltonian for one dipole excitation is given 
by 



H 



Idimer 



Ei + tefi \l,m,a) {l,m,a'\ 

a^cr' ,m,l 

+teff ^) ^ + 1> 3| + ^> 2) (/ + 1, m, 4| + H.c] . (4.6) 



Lm 



The first sum represents hopping of the dipole at a central site from one orientation to 
another, and the second sum represents hopping from one unit cell to another. Diagonalizing 



the above Hamiltonian Eq. (4.6), one obtains the spectra of a single dipole excitation given 
by 



(4.7) 
(4.8) 

(4.9) 





ky) 


= Ei- 


2tofr 




ky) 


= ^1 






ky) 


= ^1 + 


0- 




ky) 


= Ei + 





Ceff 



(4.10) 



where 2a is the distance between two neighboring central sites. The first two bands are 'fiat 
bands' as they do not depend on momentum. 
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FIG. 8: (Color online) A — )• — oo: quantum 4-state clock model: (a) decorated square lattice, (b) 
kagome lattice. Both lattices lead to the same effective Hamiltonian. Minimum of the energy for 
the effective Hamiltonian is obtained when the number of dipoles is maximum. This requires, on 
the square lattice, every central site to have a dipole, and on the kagome lattice, every site marked 
with a square to have a dipole. As a result, we can look at the system in the limit A — >• — oo as a 
collection of four state clocks, one each on the central sites of the square lattice, or on the marked 
sites of the kagome. The states on neighboring clocks are constrained by the requirement that the 
dipoles cannot overlap. 



2. Limit A — )• — oo, quantum liquid state 



To zeroth order in 1/A, all states which maximize the number of dipoles are degenerate, 
and they are all ground states of the effective Hamiltonian. These states require every central 
site on the square lattice to have a dipole (on the kagome lattice, every site marked with a 
square has a dipole), as illustrated in Fig. |8] There is an exponentially large number of such 
states which satisfy this requirement. 

To second order in 1/|A|, dipoles can hop as long as the constraint Eq. (4.2) is satisfied. 
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As there is a dipole on every central site, the constraint imphes that dipoles may only change 
their orientation while remaining at the same central site. Thus, we can write the hopping 
of the dipoles as a 4-state quantum clock model, with the clocks residing on the vertices of 
the square lattice. The hopping Hamiltonian can be written as (see Fig. Isl) 



-f^clock — Hi^ 



' l,m 
Im 

Hi^m = A - tefr(ll) + |2) + |3) + |4))((1| + (2| + (3| + (4|) (4.11) 

where (/, m) is the site index, and |cr) is a shorthand for |/, m, a) representing the orientation 
of the dipole sitting at {l,m). It should be remembered that states within the resonant 
subspace have the important constraint that two dipoles sitting at neighboring central sites 
cannot be directed toward each other. Such constraints are not contained in -ffciock- If we 
define the projection operator onto the resonant subspace which satisfies the constraints 
as P, the effective Hamiltonian with the constraints is given by projecting the hopping 
Hamiltonian, 

Hl,^^^ = PH^ioc^P. (4.12) 

In the following we will show that the exact ground state of the Hamiltonian with con- 
straint -ffciock is obtained through the projection of the ground state of the Hamiltonian with- 
out constraint, ifciock, onto the resonant subspace which satisfy the constraints Eq.(4.1). The 
ground state of the unconstraint single-site Hamiltonian Hi „i is given by |1) + |2) -|- |3) + |4) 
with eigenenergy A—At^s, and there are three degenerate excited states^ with eigenenergy A. 
Therefore, the ground state of -ffciock is given by 

iV^o) = AT JJ(|/, m, 1) + |/, m, 2) + |/, m, 3) + |/, m, 4)) (4.13) 

l,m 

where Af is some normalization constant. As we will show in the following, the exact ground 
state of H^iock i^ given by 

m=P\^o). (4.14) 

Any configuration with maximum number of dipoles corresponds to a product state given 
by Y[im\^^^^'^i,m) where ai^m = 1,2,3,4. Then the state {ipc) is nothing but the equal 
superposition of all the product states in the resonance subspace in the limit A — )• —oo i.e. 
all the configurations with maximum number of dipoles. 

The proof proceeds in two steps. First, we will show that this nodeless state \ipc) is an 



^ Interestingly, it is possible to fill the lattice with these excited states while respecting the constraint: some 
of the highest energy eigenstates of i/ciock are also eigenstates of -ff clock- Note that states of the form 



\a) — jfj'), with a ^ a', are excited eigenstates of the single site Hamiltonian (4.11). These states only 



occupy two links. We could, for example, fill the lattice with the state |1) — |2) on every site. 
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eigenstate of the Hamiltonian i^ciock- Second, we will show that it is necessarily the unique 
ground state of the Hamiltonian -ffciock- 

In order to see that \ipc) is an eigenstate of -f^ciock' check that -f^ciock l^c) is again the 
equal superposition of all the configurations with maximum number of dipoles. Consider a 
particular configuration with maximum number of dipoles represented by the product state 
|M). We count how many different configurations in produce \M) after the action of 
-^ciock- This is equivalent to counting how many different state can be created from \M) 
through the action of -f^ciock- ■^o'^ 1^)) there are in total 2N different configurations 
that are connected with \M) through an action of -f^ciock; where N is the total number 
of dipoles in the system. The decorated square lattice has in total 2N non-central sites, 
and in a configuration with maximum number of dipoles |M), of them are part of a 
dipole and of them are not. Now the action of -ffciock i^ ^ change of the orientation 
of a dipole at a central site, and the action can change the configuration by hopping a 
dipole to these non-central sites that are not occupied by a dipole in \M). Since each 
unoccupied non-central site is connected with 2 central sites, there are in total 2N different 
configurations connected with \M). This consideration is true for any |M), so we conclude 
that -ffciock iV'c) = (A — 3teff) iV'c)- Notice that the Hamiltonian Eq. (4.11) contains terms 
which bring a configuration back to itself, so the eigenenergy is (A — 3tcfr)- Indeed \ipc) 
is an eigenstate. 

Now we show that this state is necessarily the unique ground state. Note that if we take 
as the basis the product states which correspond to configurations with maximum number 
of dipoles, all the off-diagonal matrix element of -f^ciock ^^^^ negative. Moreover, repeated 
applications of -f^ciock ^^"^ connect any configuration with maximum number of dipoles to 
any other configurations with maximum number of dipoles. This statement is even true 
when the system is placed on a manifold with non-trivial topology, such as a torus. From 
these two conditions, it follows through Perron- Frobenius theoremP^ that the ground state 
of the effective Hamiltonian Eq. (4.11) is the superposition of product states with strictly 
positive amplitudes. Moreover, since two of such states cannot be orthogonal to each other, 
this ground state is unique. Now since we found the state l^/'c) which is a superposition 
of product states with strictly positive amplitudes, it follows that iV'c) must be the unique 
ground state of the Hamiltonian. 

In order to probe the properties of the ground state \ipc)^ we numerically computed 
correlations for the state in Eq. (4.14). For the correlations of local operators in a state 
which is the equal superposition of product states, the correlations can be related to the 
corresponding classical problempil. Here we study the density-density correlation (no,o^o,m) 
on the central sites, ni^rn is the density of the central site at measured relative to 

the density of the Mott insulator, so that ni^m can, in principle, take the values +1,0, — 1. 
Here +1 (—1) means that there is an extra (a missing) boson on that lattice site, and that 
this site has formed a dipole with a neighbor below or to the left (above or to the right). 
^z,m = would mean that this central site is not participating in a dipole bond; in the 
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distance (units of lattice spacing) 



FIG. 9: (Color online) Density-density correlation function {fiQ^ofiQ^rn) on the central sites (the ones 
marked by a square in Fig. |8]), obtained using transfer matrices and periodic boundary conditions. 
Correlations decay rapidly in an exponential fashion. A fit to the exponential f{m) = aexp{—bm) + 
c gives a = 0.9999, b = 2.225, c = 7.371 x 10~^, this corresponds to a correlation length of less 
than half a lattice spacing. 



limit A — 7- — oo all central sites have rii^m = =tl- We used row-to-row transfer matrices 
with rows of length up to 5 unit cells and periodic boundary conditions to compute density- 
density correlation^. Results for lattice sizes 3 x 10, 4 x 12, and 5 x 10 are plotted in 
Fig. [9j We obtain essentially the same results for different lattice sizes, and a fit of the 
correlation function to the exponential f{m) = a exp{—bm) + c gives a = 0.9999, b = 2.225, 
c = 7.371 X 10^^, corresponding to a correlation length of less than half a lattice spacing. 
Similarly, dipole-dipole correlations can also be calculated: let d^^i^m be an operator that 
projects onto states which have dipole orientation a at site (/, m). We find that correlations 
{da,o,odcT',i,m) - {da,o,o){da',i,m) decay exponentially. 

The exponential decay of the equal-time correlations suggests that the A — )■ — cxd liquid 
state has a gap to all excitations. Furthermore this liquid state is non-degenerate on a torus, 
implying the absence of topological order, and we expect the ground state to remain unique 
in the thermodynamic limit N ^ oo. Thus there is a possibility that the parent Mott 
insulator state at A — )■ oo is adiabatically connected to the liquid state at A — )■ — oo, without 
an intervening quantum phase transition. 

Ultracold atom experiments are expected to realize this quantum liquid state as long as 
the temperature is lower than the energy scale of tes oc /\E — U\. Since the magnitude of 
teff can be controlled by changing the tilting strength i?, it is likely that the temperature of 
the order of tofr is achievable in experiments. Moreover, recent experiments by Bakr et. al^ 
showed that the Mott phase contains particularly low entropy compared to the surrounding 
superfiuid, and fits with theoretical curve suggest that the Mott phase is as low temperature 
as a few nK. Therefore, together with the demonstration of lattice tiltin^^ and possibility 
to create various lattice structures through a holographic mask techniqud^, realization of 
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FIG. 10: (Color online) (a) The doubly-decorated square lattice, (b) Dense dipole configuration 
obtain in the limit A — )• — oo. Each dipole is represented by the thick line, and consists of a particle- 
hole pair. Links with two dipoles are identified as having a "dimer", and these are indicated by 
the ovals. These dimers are described by an effective quantum dimer model. 



this quantum liquid state should be possible with current technology. 



B. Doubly-decorated square lattice 

Although the quantum liquid state found in Section IV A has non-trivial entanglement, 



it does not have topological order. Here we show that by decorating the square lattice 
further, it is possible to obtain a model whose resonant subspace is labeled by the set of 
dimer packings of the undecorated square lattice. The dipole hopping terms appear in 
the perturbation theory of 1/A, and the effective Hamiltonian in this resonant subspace is 
expected to take a version of the quantum dimer modeP^. 

In this model, we consider a square lattice decorated by two sites at each link between 
lattice sites, as shown in Fig. [TO| i. Quantum spin models on a similar decorated square 
lattice were considered in Ref. [281 In the limit A — )■ — oo, the ground states in the resonant 
subspace are those with maximum number of dipoles. As before, there are a large number 
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of such states and they are all degenerate in this limit. The dipoles now come in two flavors: 
those which have one end on a central site of the square lattice, and those which reside 
exclusively on a link between central sites. Along a link between two central sites, the first 
kind of dipoles always come in pairs. We associate this pair with a "dimer" which connects 
two neighboring central sites (see Fig. [10|d). Because there can be only one dimer on each 
site for A — )■ —oo, the dimers are close-packed, and there is a one-to-one correspondence 
between the degenerate ground state manifold and the set of dimer packings of the square 
lattice. Including corrections in powers of 1/|A| we expect an effective Hamiltonian which 
acts on a Hilbert space of states labeled by the dimer coverings: this must have the form of 
the quantum dimer modeP^. In particular, a plaquette-fiip term appears at order 1/ |A|^^; 
up to this order, all diagonal terms in the basis of dipole coverings are equal and independent 
of the configurations. We leave the details of the calculations as well as possible extensions 
to a future work, and give a general discussion below. 

On the square lattice, quantum dimer models are generically expected to have ground 
states with valence bond solid order^. At critical points between different solid orders, 
topologically ordered spin-liquids can be obtainecP^'^. The quantum dimer model realized 
in this system consists of only a kinetic term, and no energy cost or gain for parallel dimers, 

i^dimer = -W ( 1 1 1 ) (= I + I =) ( 1 1 1 ) • (4-15) 

plaquettes 

This model is not critical, and so valence bond order is expected, but the precise nature 
of the square lattice symmetry breaking remains under study. An early exact diagonaliza- 
tion studj^ suggested a plaquette phase, more recent quantum Monte Carlo calculation^ 
predicted a columnar phase, and recently a mixed columnar-plaquette phase has been 
proposecP^. As A increases away from A = — oo, configurations with less number of dipoles 
become energetically favorable. These states are analogous to those in the doped quantum 
dimer modeP^'^, where other novel phases are possibld^. As the parent Mott insulator at 
A = oo is non-degenerate, it cannot be connected adiabatically to all the possible A — > — oo 
states with topological order or broken symmetry: there must be at least one intervening 
quantum phase transition. 



C. Kagome lattice tilted perpendicular to a principal lattice axis 

At last, we study the kagome lattice tilted perpendicular to a principle clXlS, clS illustrated 



in Fig. [TTj Because of the low connectivity in the lattice structure, dipole creations are 
frustrated in the limit A — > — oo and even in this limit, some sites cannot not participate in 
the creation of dipoles. As a result, there is no Ising ordering and the transverse superfluidity 
persists in the limit A — )■ — oo, as we describe below. 
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FIG. 11: (Color online) Kagome lattice tilted perpendicular to a principal lattice axis: resonant 
particle-hole excitations are illustrated. After a particle-hole pair is created, either the particle or 
the hole can hop along the direction transverse to the tilt. 



1. Limit A — )• oo 

The ground state is again the Mott insulator, and the excitations are particles and holes. 
For every particle-hole pair, either the particle or the hole is located on a transverse layer 



where it can hop, while the other is fixed and cannot move, see Fig. [TTJ The dispersion of 
the hopping quasiholes is 

ehik) = A — 2nQtcos{ka), 
while for hopping quasiparticles it is 

ep{k) = A — 2{no + l)tcos{ka). 

Therefore, the hopping quasiparticles have lower energy than the hopping quasiholes. The 
excitations are gapped for A — )■ oo, but excitations become gapless for some critical value Ac. 



2. Limit A — )• — oo 

The ground state in the resonant subspace has the maximum number of dipoles. Due to 
the geometry of the kagome lattice, some sites are frustrated in the sense that they cannot 
participate in the creations of dipoles. These sites allow the particles and holes to move along 
the horizontal direction perpendicular to the tilt, forming Luttinger liquids. As the kinetic 
energy of quasiparticles is lower than that of the quasiholes, the ground state maximizes 
the number of moving quasiparticles in a transverse layer and localizes all the quasiholes, 
as shown in Fig. [12] Each transverse layer thus realizes a "Tonks- Girardeau gas" of bosons 
at half-fillin^3H!^ a^^^j ^j^e ground state realizes a collection of uncoupled gapless Luttinger 
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FIG. 12: (Color online) A — t- — oo. Single non-degenerate ground state with half-filled Tonks gases 
along horizontal lines for A = — oo. 

liquids. These are gapless only in the limit A = — oo. Perturbation theory to second order 
in 1/A generates a dimerized hopping term, i.e. the effective hopping between two sites 
in a one-dimensional chain which share a neighbor above them becomes different from the 
hopping between sites which share a neighbor below. This breaks the translation symmetry 
along the one-dimensional layers explicitly, doubles the unit cell, and opens a gap, as the 
Tonks gases are exactly at half filling. Therefore, to order 1/A^ the one-dimensional systems 
remain decoupled, but they are no longer gapless. 



V. MODIFICATIONS FOR {Usl <^ t 

In this section we briefly discuss how the physics of tilted lattices is modified if the 



term in Eq. (1.1c) vanishes. In this case all processes (n^n) — )■ (n — l,n + 1) for general 
n are resonant at the same tilt magnitude. This greatly increases the size of space of 
states resonantly connected to the parent Mott insulator. As we will show in some cases the 
energy of this enlarged space is not bounded from below, and we then cannot use equilibrium 
methods to describe the physics. In other cases the resonant subspace remains bounded from 
below, but the phases differ from the ones found for \Us\ ^ t. We will illustrate this by 
two examples: square lattice tilted along a principle lattice direction, and diagonally tilted 
decorated square. We leave a more complete discussion for future study. 



A. Square lattice tilted along a principal lattice direction 

This case was considered in Ref. [31 However, the results there apply only if the resonant 
subspace is limited by a large {U^l. For small {U^l the physics is quite different, as we now 
describe. 
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In particular, we can no longer use equilibrium methods, as the energy of the subspace 
resonantly connected to the parent Mott insulator is not bounded from below, when A — )■ 
— oo. Once a particle-hole pair is created, the particle and hole are each free to move in 
the direction transverse to the tilt. It is then possible that two quasiparticles are adjacent 
to each other, in such a way that they can undergo a resonant transition. The system 
can resonantly reach states with an arbitrarily large number of bosons on a lattice site, as 



illustrated in Fig. 13, leading to an arbitrarily large negative energy. 
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FIG. 13: (Color online) Square lattice tilted along a lattice direction for C/3 = 0. The number in 
each square represents the number of bosons at that site. Here the occupancy of the parent Mott 
insulator is no = 1 . We illustrate resonant processes which eventually lead to the creation of lattice 
sites occupied by 4 bosons. It is possible to continue this pattern, thus we obtain any arbitrarily 
large number of bosons on a single site, through resonant processes. In the end, only a total 
number of lattice sites remain changed (shown in green). Having ns bosons on a site reduces 
potential energy nsinB — l)A/2. One could fill the lattice with such patterns. Thus the energy of 
states resonantly connected to the parent Mott insulator is not bounded from below. 



B. Diagonally tilted decorated square lattice 

This lattice remains stable: the energy of the enlarged resonant subspace is bounded 
from below. This stability is related to the lower connectivity of the lattice structure, and 
the suppression of transverse superfluidity in Section IV A Let us consider the phases in the 



limit A — > — 00. We distinguish two cases: filling factor uq = 1, and filling factor rio > 1- 
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FIG. 14: (Color online) Decorated square lattice in diagonal tilt, filling factor no = 1, C/3 = 
0. Quasiparticles are shown as dark blue circles, quasiholes as circles with yellow filling. Two 
neighboring quasiparticles can form another dipole bond (2,2) — t- (1,3) which reduces the energy 
by A. This is indicated in the figure by dotted blue lines. Quasiholes contain no bosons, and 
they cannot form another dipole bond. To minimize its energy, the system not only maximizes the 
number of dipoles on the lattice, but also maximize the number of neighboring quasiparticles. The 
set of states fulfilling these requirements maps to the dimer coverings of a square lattice, and one 
such dimer covering is shown above. 



1. Filling of parent Mott insulator: ng = 1 



Particle-hole pairs each reduce potential energy by A, and they may be arranged in 
such a way that each quasiparticle is adjacent to another quasiparticle. This allows for a 
third resonant transition, (2,2) — (1,3), which again reduces potential energy by A. As 
illustrated in Fig. 14, this always leads to a central site being occupied by three bosons. 
Two such sites are never adjacent to each other, so that an occupation number of four or 
more bosons on any site cannot be reached resonantly. The energy remains bounded from 
below. In the limit A — t- — 00 the resonantly connected subspace has a large number of 
ground states, and there is a one-to-one correspondence between this degenerate ground 
state manifold and the set of dimer packings on the square lattice. A plaquette-fiip term of 
the dimers appears in order 1/A^^. Thus we obtain a remarkable mapping of the effective 
Hamiltonian to the quantum dimer model. Unlike the large f/3 case in Section IV, here the 
dimer model appears already for the singly-decorated square lattice, and so should be easier 
to realize experimentally. 
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FIG. 15: (Color online) Decorated square lattice in diagonal tilt, filling factor tlq > 1, U3 = 0. In 
this case 'quasiholes' can also undergo a resonant transition, which is shown as orange dotted lines. 
There are two degenerate ground states with density-wave order: the one shown in the figure, and 
its symmetry related partner. 



2. Filling of parent Mott insulator: uq > 1 



If the filling factor is two or larger, then two adjacent 'quasiholes' can undergo a resonant 
transition (no — 1, no — 1) — )• (no — 2, no), reducing the potential energy by A. Potential 
energy is minimized if each 'quasihole' and each 'quasiparticle' can undergo such a transition, 
as illustrated in Fig. 15 This leads to a two fold degenerate ground state with density- wave 
order: all central sites are occupied by either (no + 2) or (no — 2) bosons, while all other 
sites contain no bosons. 



VI. CONCLUSIONS 



This paper has shown that there are rich possibilities for generating non-trivial quantum 
states upon tilting a Mott insulator of bosons. Our classical intuition tells us that applying 
a strong tilt to a Mott insulator should lead to a runaway instability of particles flowing 
downhill. However, quantum mechanically, for a single band model in the absence of external 
dissipation, and with a significant three-body interaction U^, this does not happen. Instead 
particle motion is localized in the direction of the tilt, and free motion is only possible in 
orthogonal directions. It is useful to introduce the idea of a resonant subspace of states 
which are strongly coupled to the parent Mott insulator when the tilt approaches its critical 
value. Within this resonant subspace, we can define an effective Hamiltonian whose energy 
is bounded from below, and which has well defined equilibrium quantum phases and phase 
transitions. We presented a variety of such effective models here, and described general 
aspects of their phase diagrams. As in the previous worlP, we found phases with Ising 
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density wave order along the tilt direction, and superfiuid/Luttinger liquid behavior in the 
direction transverse to the tilt. More interestingly, we showed that on a variety of frustrated 
or decorated lattices, particle motion can be suppressed also in directions transverse to the 
tilt. Then, the tilted lattices map onto quantum clock or dimer models, with novel quantum 
liquid and solid phases. 

We also briefly noted, in Section |V| the situation when \U3\ was small. In some cases, the 
energy of the resonant subspace is unbounded from below, and so non-equilibrium methods 
will be necessary to understand the physics. However, for the simple case of the decorated 
square lattice, we found a mapping of the effective Hamiltonian to that of the quantum 



dimer model in a limiting illustrated in Fig. Q4l 

It is clear that the models studied here have extensions to numerous other lattices and 
decorations, and that there are many promising avenues for obtaining exotic phases. It 
would also be interesting to study the analogous models of fermions, which could display 
metallic states associated with motion in the transverse directions. In principle, all these 
models should be readily accessible in experiments on trapped ultracold atoms. 
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Appendix A: Decorated square lattice in the A — t- — oo limit 
1. Ground state correlation functions 

To calculate correlation functions in the exact ground state of the effective clock model, 



Eq. (4.14), we make use of the following "transfer matrix" method. For an equal-amplitude 
ground state, correlations of operators, O, which are diagonal in the basis of dipole coverings, 
can be related to a corresponding classical problem 

,^,_ {^Pc\dm _ j:M{M\d\M) 



26 



Here, the sum over M runs over all states that maximize the number of dipoles and respect 



constraints Eq. (4.1) and (4.2). We assume that all the states |M) are properly normalized, 
(M|M) = 1, and define = J2m(^\^) = dim(ifc) as the dimension of the Hilbert 
space with constraint. The density operator ri; ^ measuring the relative density at central 
site {l,m) is diagonal in the basis of dipole coverings; it has eigenvalues ±1. Therefore, 

(^(o,o)^(/,m)) can be written as 

(^(o,o)^a,m)) = , (A2) 

where N± is the number of configurations which have (M| n(o,o)^{;,m) \M) = ±1. The problem 
of calculating correlation functions of the ground state reduces to counting the number of 
classically allowed dipole coverings. This counting can be done using row-to-row transfer 
matrices. 

Let T be the transfer matrix for a row of length A^^ unit cells. T is a 4^"= x 4^"= matrix 
and acts on the space of dipole configurations within that row. Note that this Hilbert space 
includes configurations which are not allowed by the constraint. Such configurations are 
excluded by setting the corresponding matrix elements of T to zero. The matrix elements of 
T are either 1 or 0, so that T^/ c = 1 if the sequence of two row configurations c, c' contains 
no violations of the constraint, and Tc'^c = otherwise. The dimension of the Hilbert space 
for a Nx X Ny lattice with periodic boundary conditions is then given by the trace of powers 
of the transfer matrix, 

= Tr(T^^). (A3) 



A naive numerical implementation of Eq. (A3) (without using the sparseness of T) gives a 



complexity of 0(A^y4^^"=), compared to going over all the states in the Hilbert space with 
complexity 0(4^^^^). The values for N± can be similarly found by including projection 
matrices in the trace, for example 

A^+ = Tr {T^y~'^Pp.iT"'Pp.o) + Tr (T^^-^P^.^T^P/^.q) . (A4) 

Here Pp-i {Ph-i) is a projection matrix that projects onto configurations which have a particle 
(a hole) at site / = 0, 1, . . . , A^^^ — 1 of a given row, respectively. Using this method different 
correlations for the dipole-directions can also be calculated. 



2. Connecting the ground state to a site-factorizable state 



In this appendix, we provide evidence that the ground state of -ffciock) -^l- (4-14), is not 
topologically ordered. To show this, we demonstrate that the ground state of -f^ciock "^^^ 
be smoothly connected to the ground state of -ffciock (without the "hard core" constraint of 
the dipole-clocks). The ground state of -f^ciock is a site factorizable state, and therefore it is 
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topologically trivial. The fact that l^/^c) can be connected to it, while keeping the correlation 
length finite, is a strong indication that \ip(^ does not have topological order either. To 
connect the ground state of -ffciock ^^e ground state of i^ciock, Eq- (4.11), we define the 
wave function 

|^(e)) = 5^6"(^)/^|C). (A5) 

c 

Here the sum over C runs over all 4^ dipole coverings with one dipole per unit cell, including 
those which do not respect the constraint. The function n{e) is an integer counting the 
number of colliding arrows in configuration C. All allowed configurations M have n{M) = 0. 
Thus changing the value of e interpolates between the ground state of -f^ciock' l^c), and the 
ground state of i^ciock, iV'o), 

\^{e = 0)) = IV^,) , (A6) 
me = 1)) = m) . (A7) 

As the density operator „j) is still diagonal in the basis of C, we can again use (generalized) 
transfer matrices to compute density-density correlation functions for any value of e. The 



result is shown in Fig. 16 , and we find no divergence in the correlation length. This supports 
our conjecture that the ground state \ipc) is topologically trivial. Of course, we have not 
calculated all possible correlation functions, so this is an evidence, but not a proof. 



To calculate the correlation functions, we make use of equations (A2), (A3), and (A4), 
with one modification: the transfer matrix T(e) now depends on the parameter e in the 
following way 

(T(e)),,e = (A8) 

where n(c', c) is the number of collisions in row-configuration c', plus the number of collisions 
between rows c and c' . The number of collisions in the other row, c, does not come in here, 
to avoid double counting. 
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Interpolating density-density correlations 




distance in units of lattice spacing 



FIG. 16: (Color online) Density-density correlations of a wavefunction which interpolates between 
the ground states of H^^^^^ and Hdock, see text. Different colors stand for different values of the 
parameter e. The correlation length does not diverge, it decreases as e — )• 0. This provides evidence 
that the state \ipc) is topologically trivial. 
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